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We derive a simple closed analytical expression for the total entropy production along a single 
stochastic trajectory of a Brownian particle diffusing on a periodic potential under an external 
constant force. By numerical simulations we compute the probability distribution functions of 
the entropy and satisfactorily test many of the predictions based on Seifert's integral fluctuation 
theorem. The results presented for this simple model clearly illustrate the practical features and 
implications derived from such a result of nonequilibrium statistical mechanics. 
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A great deal of interest has been recently devoted 
to the development and understanding of formal rela- 
tionships known as fluctuations theorems. These ex- 
pressions generally constrain the probability distribution 
function of entropy production without relying on how 
close the system is to thermal equilibrium. A num- 
ber of complementary proposals have been put forward 
0, B 1 S IS S, 0, S, (1 E3, and a considerable effort 
has been devoted to clarify the similarities among them. 
Experiments have been carried out to validate and use 
such expressions to extract useful information in differ- 
ent physical and biological systems [II], [l2[ • A new fluc- 
tuation theorem for dissipative systems has been derived 
by Seifert [13], which is based on the generalization of 
the concept of entropy production to the level of a sin- 
gle trajectory. The overall entropy production, namely 
the sum of the contributions from the environment and 
the one from the particle along its trajectory, is shown to 
obey an integral fluctuation relation (and in some cases a 
detailed one) for any initial conditions, arbitrary driving 
and over a generic finite time interval. Moreover, recent 
experiments have validated this theorem [3, [IB] and the 
stochastic entropy production has been measured exper- 
imentally for the first time [Tfl ]. 

In this work, we consider a Brownian particle moving 
on a periodic potential under the action of an external 
constant force. Such system allows us to carry out a de- 
tailed study of the total entropy production and assess 
carefully the implications of Seifert's fluctuation theorem. 
After a summary of the main results of Ref. [IH in Sec- 
tion HIl we particularize them for the system under study 
in Section IIIII In Section IIV1 we discuss numerical re- 
sults and compare them with analytic predictions where 
available. We conclude with a discussion section. 



In Ref. [l3[ the dynamics of one degree of freedom x 
such as the position of a particle in a one dimensional 
system is studied. It evolves according to a Langevin 
equation, 



x = nF(x,\(t)) +X; 



(1) 



where x is a Gaussian white noise with autocorrelation 
(x(t)x{t')) — 2D5(t — £'), and D is the diffusion constant 
so that in equilibrium D = T[i. We consider fog = 1 
throughout this work. The force F(x,X(t)) = — V'(x) + 
4>{x, X(t)) splits into a term arising from a conservative 
potential, V(x), plus a time dependent external source, 
4>(x,X(t)), whose control parameter is X(t). The non- 
equilibrium Gibbs entropy, 



S(t) 



dxp{x,t)\np(x,t) = (s(t)}, 



(2) 



can be interpreted as a mean over single trajectories as 
specified by the probability p(x, t) of finding the particle 
at position x at time t. Seifert uses the above expression 
to propose an entropy associated to each trajectory of 
the particle, 



s(t) = -lnp(x(t),t). 



(3) 



The rate of change of entropy of the system is de- 
rived from the Fokker-Planck equation for the probability 
p{x,t). The rate of heat dissipation in the medium and 
the increase of entropy in the medium can be related as 
q(t) = F(x, X)x = Ts m . Therefore, it is possible to write 
down the rate of change of total entropy associated with 
a single trajectory or realization as 



Stot{t) = s m (t) + s(t) = - 



d t p(x, t) 



p(x,t) 



X (t) D P(x,t) 



x(t) 



where we have introduced the particle flux, which for the 
Langevin dynamics we are considering is given by 
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In the spirit of previous fluctuation theorems, which 
relate the weights of trajectories with their time reversed 
counterparts, Seifert introduces the quantity R, 



R[x(t),\(t);po,pi] 



In 



p[x(t)\x ]p (x ) 
p[x(t)\x ]pi(x ) ' 



(6) 



which relates the weight of a trajectory with its time 
reversed analogue. The (arbitrary) distributions of ini- 
tial and final values are Po(xq) and Pi(xq) respectively. 
The tilde refers to the time reversal operation, which as- 
sociates with a trajectory x(t), a reversed one x(t) = 
x(t — t), where t is the duration of the process. In this 
way, x(t) retraces the reversed path of the stochastic tra- 
jectory x(t), and it starts at Xq = 5(0) = x{t), which 
corresponds to the final value of the stochastic trajec- 
tory x(t). The system is driven out of equilibrium by an 
external control parameter that changes with time ac- 
cording to a prescribed protocol X(t), and p[x(£)|x ] and 
p[x(i)|xo] are the probabilities of the forward and back- 
ward paths respectively. 

While R is in principle defined for any normalized dis- 
tribution pi(xo), it only becomes the total entropy pro- 
duction (and thus acquires a physical meaning) when 
Pi(x t ) =p(x,t) is the solution of the Fokker-Planck equa- 
tion for the initial distribution p (xo). Then, the defini- 
tion © yields As = lnp (xo)/pi(x t ), the particle en- 
tropy change. The other terms appearing in Eq. ([6]), re- 
lated to the path probabilities, correspond to the medium 
entropy production As m = lnp[x(i)|xo]/p[x(t)|xo], as 
shown in Ref. 17]. Finally, we can identify R with the 
total entropy production of the system, 



R = As + As r 



Asto 



(7) 



As we will show in short, the mathematical identity 
(e~ R ) = 1, together with Eq. ([7]), implies the following 
integral fluctuation theorem for the total entropy varia- 
tion along a stochastic trajectory Astot, 



-As to 



1. 



(8) 



It is absolutely general since it holds at any time along a 
trajectory generated from generic initial conditions and 
subject to arbitrary forces. In the subsequent sections 
we will analyze the implications and consequences of this 
relation for a simple example where a detailed analysis 
of R will be performed. 



III. DERIVATION OF AN EXPRESSION FOR R 

Let us consider a one dimensional stochastic system in 
which a particle is moving in a periodic potential V(x) 
and under an external force F (so that we restrict the 
general driving to the simplest case of constant driving 
in time and space, <f>(\(t),x) — F) . The corresponding 
Langevin equation in the overdamped limit is 



x = -V'(x)+F + ri(t), 



(9) 



where the friction coefficient has been absorbed in the 
time units for simplicity. Then the noise is Gaussian, 
white, with zero mean and autocorrelation (rj(t)ri{t')) — 
2TS(t — t'). This system is a simplified version of mod- 
els analyzed previously in Refs. [H, [l9[ to assess some 
aspects of the fluctuation theorem. 

Our aim is to derive an explicit expression for R in 
the steady state. To this end, we first split R into two 
contributions, 



R=ln p\mxo\ +ln po(*o) 



p[x(t)\x ] 



(10) 



The first term on the rhs involves the weight of the for- 
ward and reversed paths. From the path integral ap- 
proach to stochastic systems, it is known that the prob- 
ability to observe a given realization of the noise is 



p[r/(t)] ~ exp 



dt'r) 2 (t')/4T 



(11) 



Such a Gaussian dependence can be used, together with 
the evolution equation (Eq. (9])), to obtain an explicit ex- 
pression for the transition probability p[x(i)|xo] in terms 
of the particle's position. If one takes into account the 
Jacobian of the transformation, one arrives at [l7| 



In 



p[x(t)\x ] 



(-V'(x)+F)xdT/T, (12) 



' p[x(t)\x ] 

which in the present case reduces to a much simpler form, 



In 



P[x(t)\x ] 
p[x(t)\x ] 



(FAx — AV)/T, 



(13) 



where Ax = x(t) — x(0) = Xt — xo and AV = V(x t ) — 
V(xq). It is worth mentioning that the final position x t 
corresponds to a generic observation time t. 

The second term on the rhs of Eq. (fTU|) accounts for 
the ratio of the probabilities of initial and final positions. 
We are interested in the steady state, and thus p{x) fol- 
lows easily from the corresponding solution of the Fokker- 
Planck equation. For a system under an effective poten- 
tial U(x) = V(x) — xF, being V(x) an arbitrary periodic 
function of period L, one gets [2(| 



p(x) = {M/T)e- U( - X)/T 



dye 



U(y)/T 



(14) 



where Af ensures proper normalization. Since the system 
is in the steady state, po(x) = Pi(x), which implies 



Pojxo) _ p{x Q ) 
Pi(x ) p(x t )' 



(15) 



and hence the second term of Eq. (jT0|) can be obtained 
from Eq. (Til]) , 

When inserting Eqs. (13]), and $U$) plus (15]) into (fTQf , 
a simple expression for R is obtained, 



rX t +L , 

L t d v e 



R = In 



U(y)/T 



Q +L dye u (v)/ T 



(16) 
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Figure 1: A) Periodic potential defined in Eq. (|17[1 for L = 2-7T. 
B) Total entropy variation R obtained from Eq. (|19|l as a 
function of the final position xt (when the initial position is 
xo = 0), for different values of the dimensionless external force 
parameter /. 



It shows that R can be computed exactly by integrating 
the effective potential felt by the particle for every trajec- 
tory. It is worth mentioning that this expression is only 
a function of the initial and final positions, which is a 
non trivial result. This fact is of experimental relevance, 
since one does not need to keep detailed information of 
the complete trajectory to compute R. Note also that 
the exchange of xq for xt reverses the sign of R. We will 
exploit this important symmetry subsequently. 

Further progress is feasible if a specific potential is cho- 
sen, so that the integrals determining R in the previous 
equation can be worked out analytically. In particular, 
by picking the potential depicted in Fig. (JTJ A), 



V(x) = V ln(cos(27ra;/L) + 2), 
Eq. (fl~6)) reduces to 



R = ln Jxt 



dye-yi (cos(y) + 2) 1 



Q +2 "dye-yf (cos( 2/ ) + 2) 1 



(17) 



(18) 



where / = FL/2irT and v = Vq/T. For natural values of 
the parameter vq, the integrals can be done analytically, 



yielding to an explicit expression for R as a function of 
the initial and final values of the position of the particle 
for a trajectory taking place during an interval t. There- 
fore, and without loss of generality, we consider vq = 1 . 
We then arrive at the simple expression 



R = f(x t - x ) - In 



Kxt) 



where 



I( x ) = f 2 (cos(x) + 2) - /sin(x) + 2. 



(19) 



(20) 



The first term on the rhs of Eq. (fl9|) is directly related 
to the external work performed on the system by the ex- 
ternal force F and then dissipated. The second term has 
information about the potential V(x) and of the external 
force too. In Fig. ([IjB) we show the increase of R as a 
function of the final position along a trajectory, xt, when 
xq = for different values of the dimensionless driving 
force /. When / is increased, the relative contribution of 
the periodic potential to R decreases, as can be seen di- 
rectly from Eq. (|19p . The interesting regime corresponds 
to the situations in which the external driving, the peri- 
odic potential and the thermal energy have comparable 
magnitudes. We will concentrate on such regime from 
now on. 



IV. NUMERICAL RESULTS 

We perform numerical simulations of Eq. @ with the 
potential defined in Eq. (|17p to characterize the proper- 
ties of R. Such study is required because, even though R 
depends only on the initial and final positions of the tra- 
jectory, the final position is a function of the noise and 
hence the probability distribution function of R 7 p{R), 
has to be obtained. Then its properties can be analyzed 
and the fluctuation theorem tested. 



A. Integral fluctuation theorem 

From the definition of R one can show that 



x(t),X 



p[x(t)\x ]p (x )e = 



^2 P[£(*)l*o]Pi(£o) = !• 

i(t),x 



(21) 



This is a general result, independent of the specific form 
of the potential and external driving exerted on the par- 
ticle. Using the collection of values of xq and Xt from the 
simulations and inserting them in expression (|19p , we can 
perform the average (e~ R ), which can be written down 
as the sum 



1 N 

N ^ 

i=i 



(22) 
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Figure 2: Average of e~ fl as the number of trajectories N in 
the statistics is increased. The parameters are / = 1 and t = 
1. Convergence is reached quickly. Sudden leaps are observed 
as shown in the inset, which is a zoom of the main figure. 
Such jumps correspond to the occurrence of rare realizations 
(R < 0), which are crucial contributions to the average. 



where N is the number of particle trajectories that are 
used to perform the average. It can also be interpreted 
as the number of repetitions of the experiment with the 
same protocol. We display in Fig. @ the average of 
e~ R as a function of the number of trajectories, when 
the time interval is t = 1 and the external driving is 
/ = 1. Around this zone of the space of parameters is 
where the relation is not trivial at the same time that it is 
numerically testable. One can see that the average value 
predicted by Eq. (PT|) is rapidly attained. The same must 
hold for any values of vq, f and t, as it has been verified 
for some of them (data not shown) . Although a noisy but 
progressive approach to unity is found as N increases, 
jumps are observed, as depicted in the zoom of Fig. $Z§. 
They are rare events that contribute significantly to the 
average because R < 0, and hence can be regarded as 
transient violations of the second law. We will come back 
to this point in Fig. ^ when the evolution of R in time 
for one single trajectory will be discussed. 



B. Probability distribution function of R 

In order to gain more insight into the properties of R 
and the implications behind the fluctuation theorem, we 
analyze the probability distribution of the total entropy 
production, p(R), as obtained from numerical simula- 
tions over many independent trajectories. Such a quan- 
tity can be measured experimentally, and a specific real- 
ization of a driven two level system with time-dependent 
rates has been already achieved [l4T |. 



Figure 3: Probability distribution functions of R at different 
times t for an external driving / = 1. 



Fig. ([3]) exhibits the probability distribution function 
p(R) at different times, for an external driving / = 1. 
The distribution function is peaked and would reduce to 
a delta function at R = at all times for an equilibrium 
system. Due to the increase of the mean of R as a func- 
tion of time, p(R) is shifted towards positive values of 
R and spreads as time increases. The broadening of the 
distribution can be attributed to the interaction of the 
particle with the underlying periodic potential together 
with thermal fluctuations. It is interesting to note that 
p(R) develops a richer structure as time increases. In 
particular the number of peaks increases because there 
is a higher chance to explore the periodic potential. In- 
creasing / would lead to a corresponding decrease in this 
structure in p(R). For short time intervals, there is a 
contribution to the distribution that comes from nega- 
tive values of R, in agreement with the jumps described 
in Fig. |(2J). As time increases the contributions to the 
negative tail decrease, being harder to observe. Precisely, 
the relation between negative and positive tails is very 
relevant, as we will discuss in short. 



C. Decomposition of R 

Let's recall that R can be expressed as a sum of two 
contributions, as written down in Eq. (| 10[) . The first term 
involves a ratio of conditional probabilities and is inex- 
tricably related to the heat Q interchanged with the heat 
bath (p, because Q/T = As m = f*W(x-r)(t'))dx(t')/T 

reduces to Eq. (fl"3|) when Eq. ([9]) is used. In turn, Q can 
be decomposed in two significant contributions; a first 
term proportional to the work produced by the external 
force, W = FAx/T, and a second one, V = -AV/T, 
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Figure 4: Probability distribution functions (at t = 5 and 
/ = 1) of the several contributions to the total entropy change 
R: A) Total entropy production R (solid line), dissipated 
heat Q (dashed) and work W (dotted). B) Internal energy V 
(dashed) and particle entropy change As (solid). Such quan- 
tities satisfy R = Q + As and Q = W + V. 



Figure 5: p(R) for / = 1 at t = 10. The dotted area shows the 
part of the probability distribution function used in the inset. 
The second peak is not included in the range of the inset 
because the statistics of the corresponding negative values 
is very poor. The detailed relation \np(R)/p(—R) = R is 
checked satisfactorily. 



which corresponds to the change in the internal poten- 
tial energy of the system. The remaining contribution 
to R, the second term on the rhs of Eq. (|10p . is the 
entropy change of the particle As along the stochastic 
trajectory. On physical grounds, it is fruitful to write 
R = Q + As = W + V + As so that three different 
physical sources contribute to the change in the overall 
entropy of the system. The total entropy change R has 
contributions from the medium and from the system it- 
self, namely, the total heat dissipated Q (which comes 
both from the external forcing and the change of internal 
energy) and the entropy change associated to the change 
in the system's internal state, As. 

As presented in Fig. ([5]), for the simple model under 
scrutiny it is possible to study each probability distri- 
bution function of the different contributions separately. 
Obviously, although R = Q + As, p(R) ^ p(Q)+p(As). 
One can see that the structure of the probability dis- 
tribution function p(W), whose contribution is directly 
the external force, has no internal structure and is sim- 
ply characterized by a broad, slightly tilted, peak which 
moves to higher values of W as time increases. p(Q) re- 
tains the structure of the peaks, which we attribute to the 
effect of the periodic potential. See Ref. [22| for a study 
of the structure of the probability distribution function of 
heat exchanged in a simple temperature transducer. In 
Fig. ([UB), p(V) is shown. It has a bound support and its 
shape does not change as time increases. Since the term 
associated to the external force takes increasing values, it 
is clear that the contributions from the periodic potential 
to the entropy change will have a restricted effect on the 



overall p(R). Nevertheless, the potential is the responsi- 
ble for the peaked structure in p(R). The contribution 
of As is qualitatively analogous to that of V. 



D. Detailed fluctuation theorem 

A stronger fluctuation relation has also been shown to 
hold for the probability distribution function of R itself 

EE!, 



In 



Kg) 

p{-R) 



= R. 



(23) 



It constrains directly the properties of the probability dis- 
tribution function of R. This stringent relation is cl osely 
related to the studies performed in Refs. [H, 0, 0, llOj . 
where similar fluctuation relations are derived for en- 
tropy, work and heat in different contexts. The previous 
relation follows [23| from the fact that for an arbitrary 
function g(R), 

(g(R)e- R ) = ]T p[x(t)\x ]p (x ) g(R)e- R = 

x(t),x 

= J2_pi S; ( t )\^}pi(xo)g(-R[x(t)}) = {9(-R)), (24) 

x(t),x 

where the symmetric property of R, namely i?[a;(t)] = 
—R[x(t)], is used and necessary to prove the above equa- 
tion. This important condition must be verified in every 
case. Here it is satisfied because R reverses sign when 
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Figure 6: Test of the relation ([25} for the zeroth, first, second 
and third moments of p(R) at t = 1 and / = 1 as a function 
of the number of trajectories. 



interchanging cc an d x t . Then, for the particular case 
g{R) = 5(R[x(t)\ - R), one recovers Eq. Ip3jl . 

In the inset of Fig. ([5]) we check the detailed fluctua- 
tion theorem by relating the negative and positive con- 
tributions to the probability distribution of R, leading 
to an excellent agreement with the theoretical predic- 
tion. As we will show in detail, (R) increases linearly 
with time. Therefore, the probability distribution func- 
tion of R moves towards positive values as time goes by. 
Then, concerning the verification of the detailed fluctu- 
ation theorem, the longer the time of observation, the 
poorer the statistics at the negative tails, thus limiting 
the possibility of a meaningful tests. Experimental tests 
have been performed in Ref. [lj| to verify the detailed 
fluctuation relation. 

One can further take profit of Eq. by considering 
g(R) = R k . A set of relations involving average moments 
of powers of the total entropy change of the system are 
found, 



(R k 



(-l) k (R k ) 



(25) 



In Fig. ([6]) we display the lower order moments (k — 
0,1,2, 3) as a function of the number of independent tra- 
jectories, N. Theoretical predictions are recovered. The 
dependence on N shows again the statistical nature of 
the fluctuation theorem. Note that the behavior is anal- 
ogous to the one in Fig. 



Figure 7: Probability distribution function p(R) for / = 1 at 
t = 1 obtained from different ensemble initial conditions: cir- 
cles correspond to steady state initial distribution of positions, 
triangles to a delta function at position zero and squares are 
obtained from a Gaussian distribution of positions. The in- 
set shows that only the probability distribution function of R 
obtained from the steady state ensemble of initial conditions 
fulfills the detailed version of the theorem, Eq. (|23|l . 



only in the steady state. Therefore, when sampling from 
other different initial conditions, the relations verified so 
far will not be reproduced anymore. Fig. ([7]), displays 
the histograms at t = 1 corresponding to the dynamics 
of the system when three different ensembles of initial 
conditions for the position are considered: steady state, 
Gaussian distributed and a delta function at zero po- 
sition. The inset shows the ratio lnp(R)/p(— R) as a 
function of R for all three initial distributions. One can 
clearly see that only for the steady state initial ensemble 
the detailed fluctuation theorem holds. It is worthy to 
stress this point to prevent confusions and errors when 
applying the predictions provided by the theorem in real 
experiments, in which one should have careful control on 
the preparation of the initial state of the system. It would 
be interesting, although more tedious, to test the predic- 
tions derived from the fluctuation theorem away from 
the steady state, in which no simple expression such as 
Eq. (TIT)]) is available. 



F. Mean dissipated work and (R) 



E. Initial conditions 

Seifert's integral fluctuation theorem is valid for arbi- 
trary initial conditions. However, Eq. (fT!?)) holds for R 



For dissipative systems in a steady state, R is closely 
related to the so called housekeeping heat Qhk . This is 
the heat that is permanently dissipated while maintain- 
ing a nonequilibrium state at fixed external parameters. 
It has been shown [l8| that R — Qhk/T, at the same 
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time that 



(Q 



dx 



(26) 



The steady flux J is related to the probability p(x) ac- 
cording to Eq. ([5]) , which makes it possible to express the 
housekeeping heat in terms of the external driving and 
the particle flux as (Qhk) — tFJL. Since JL = (x), the 
mean heat produced can be rewritten in the suggestive 
form 



(Qhk) = F(x)t, 



(27) 



which implies that the mean total entropy change can 
be simply expressed in terms of the mean work produced 
and dissipated by the driving force, (W) = F(Ax), as 



(R) = (W)/T, 



(28) 



We can derive an explicit expression for the mean velocity 
of the particle [2(| 



LT{\ - e- 2 *f) 



J L dxe- u ^)' T f* +L dye u (v)/ T ' 



(29) 



Remember that U(x) = V(x) — Fx. Again, for the po- 
tential V(x) given by Eq. Q17p with vq = 1, one gets, 
after some integrations, 



2ttT p + f 
L P + 2/V3' 



(30) 



Note that the mean velocity could be recasted in terms 
of an effective mobility (i e f f , which would depend on the 
applied force, so that (x) — F/i e ff, For small forces, 
/i e // = y/3/2 which would correspond to the free diffu- 
sion of a Brownian particle in the periodic potential, and 
for asymptotically large forces, fJ- e ff reduces to the bare 
mobility, fi e ff = 1, in scaled units. 

Finally, coming back to the mean entropy production, 
the following analytical expression is then encountered, 



(R) = tuf 



/ 2 + l 

/ 2 + 2/V3' 



(31) 



where u> = T(2ir/L) 2 . In Fig. © we compare the theo- 
retical prediction for (R) with numerical results, explor- 
ing the force / and the time t. Theory and simulations 
are in excellent agreement. The linear dependence in 
time of the mean value of R yields to a constant mean 
entropy production rate. 




Figure 8: Plot of (R) versus the external force / for tu) = 1 
(A) and versus the dimensionless time tu for / = 5 (B). Dots 
are obtained from numerical simulations, while solid lines cor- 
respond to the exact analytical expression (|31l> . 



thermal energy, negative values of R will be observable, 
corresponding to the so-called "transient violations of the 
second law" . These violations contribute significantly to 
the averages at the same time that determine the tails of 
the distributions, which are crucial when examining the 
detailed fluctuation relation. In fact, negative values of 
R are directly related to Figs. ((2]), ([3|) and |5]), and their 
corresponding subsections. 

In Fig. @ we plot the time evolution of the entropy 
production R of a single trajectory. One can see that 
at short times R is likely to become negative, indicat- 
ing that up to that point the total entropy production 
associated to that stochastic trajectory is negative. The 
inset shows that at long times R is positive and increases. 
Upon averaging over a number of independent trajecto- 
ries, the linear increase predicted in the previous section 
is recovered. Hence, the study of single trajectories high- 
lights the stochastic nature of fluctuation theorems and 
the need to distinguish between mean response of the 
system and the properties of individual trajectories. As 
the period of observation increases, it becomes rapidly 
unlikely that negative values of R are observed. This, 
together with the energy scale, sets the regime where ex- 
amining "transient violations" becomes impractical. 



CONCLUSIONS AND FINAL REMARKS 



G. "Transient-violation" trajectories 

Although (R) increases linearly in time, for a single 
trajectory R will fluctuate. As a result, for short times 
and for typical energies of the system comparable to the 



We have considered a simple system in a nonequilib- 
rium steady state and have studied in detail many of the 
predictions and consequences derived from the fluctua- 
tion theorem in Ref. [13j |. In particular, by choosing an 
appropriate periodic potential, we have been able to find 
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Figure 9: Total entropy change along one single trajectory Ri 
as a function of time for f — 1. Note that, for short times, 
"transient violations of the second law" are observable since 
the total R is negative. In the inset the behavior of R at 
larger times is plotted, where fluctuations are less dominant. 
Averaging over a large number of independent trajectories 
would lead to the linear increase of (R) predicted by Eq. (|3T)l . 



an explicit analytical and physically meaningful expres- 
sion for the total entropy production along a stochas- 



tic trajectory. Advantage has been taken from the fact 
that only initial and final positions for each trajectory 
are relevant, being a potential useful feature in future 
experiments. Numerical simulations of the dynamics of 
the Langevin equation have been necessary to build the 
probability distribution functions of R. We have eluci- 
dated the relevance and properties of p(R) in the system. 
Moreover, all of the results based, related, derived and 
predicted from the fluctuation relation have been tested 
satisfactorily. 

This work leads to a more profound understanding of 
the specific meaning and practical appliance of such ab- 
stract general results. At the same time, it opens the 
door to study new and more complex cases such as inter- 
acting particles or the effect of temperature dependent 
substrates. Furthermore, it would be very interesting to 
analyze fluctuations observed in biological systems such 
as molecular motors under the perspective of fluctuation 
relations. Studies based on simple physical models are 
crucial to illustrate the key ingredients and foresee the 
consequences of these theoretical results on real experi- 
ments. Therefore, they bridge the fundamental advances 
in statistical mechanics with the active area of nano scale 
experiments on biological systems. 
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